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We present a theoretical study of the Goos-Hanchen and Imbert-Fedorov shifts for a fundamental Gaussian 
beam impinging on a surface coated with a single layer of graphene. We show that the graphene surface 
conductibility rr(uj) is responsible for the appearance of a giant and negative spatial Goos-Hanchen shift. © 
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When an optical beam impinges upon a surface, non- 
specular reflection phenomena may occur, such as the 
Goos-Hanchen (GH) [1-3] and Imbert-Fedorov (IF) [4,5] 
shifts, resulting in an effective beam shift at the interface. 
A comprehensive review on beam shift phenomena can 
be found in Ref. [6]. Although Goos and Hanchen pub¬ 
lished their work more than 60 years ago [1], this field 
of research is still very active, and in the last decades 
a vast amount of literature has been produced on the 
subject, resulting not only in a better understanding of 
the underlying physical principles [7 -14], but also in the 
careful investigation of the effects of various held config¬ 
urations [?, 15,17-19] and reflecting surfaces [20-23] on 
the GH and IF shifts. 

In recent years, on a parallel trail, graphene attracted 
very rapidly a lot of interest, thanks to its intriguing 
properties [24, 25]. Its peculiar band structure and the 
existence of the so-called Dirac cones [26], for exam¬ 
ple, give the possibility to use graphene as a model to 
observe QED-like effects such as Klein tunneling [27], 
Zitterbewegung [28], the anomalous quantum Hall ef¬ 
fect [29] and the appearance of a minimal conductiv¬ 
ity that approaches the quantum limit e 2 /fi. for vanish¬ 
ing charge density [30]. In addition, the reflectance and 
transmittance of graphene are determined by the fine 
structure constant [31], and a single layer of graphene 
shows universal absorbance in the spectral range from 
near-infrared to the visible part of the spectrum [32]. 

Among the vast plethora of applications, graphene also 
proved to be a very interesting system where to observe 
beam shifts. Very recently, in fact, the occurrence of GH 
shift in graphene-based structures has been reported, 
both for light beams (where giant GH shift has been 
observed [33]) and for Dirac fermions [34]. Despite all 
this, a full theoretical analysis of GH and IF shifts in a 
graphene-based structure has not been yet carried out. 

In this Letter, we therefore present a theoretical anal¬ 
ysis of the GH and IF shifts occurring for a monochro¬ 
matic Gaussian beam impinging onto a glass surface 
coated with a single layer of graphene. The results of our 
investigations show on one hand, that the appearance of 


a giant GH shift is ultimately due to the graphene’s sur¬ 
face conductivity cr(w), and on the other hand, that the 
presence of the single layer of graphene introduces a de¬ 
pendence of the phases of the reflection coefficients on 
the incidence angle, thus resulting in a nonzero spatial 
GH shift also when total internal reflection does not oc¬ 
cur. 

We start our analysis by considering a monochromatic 
Gaussian beam with frequency ui = ck (with k being the 
vacuum wave number), impinging on a dielectric surface 
characterized by the refractive index n and coated with 
a single layer of graphene [Fig. 1 (a)]. The graphene layer 
is characterized by the optical conductivity cr(tu), whose 
expression can be given in the following dimensionless 
form [26] 


Act 

<J(U) — i— + 7ra 


e(n - 2) + ^ ln PTf' 


, ( 1 ) 


where Q = is the dimensionless frequency, y, is 

the chemical potential, a ss 1/137 is the fine structure 
constant [28] and 0(a:) is the Heaviside step function 
[35]. 

According to Fig. 1(b), we define three Cartesian ref¬ 
erence frames: the laboratory frame K = ( 0,x,y,z ) at¬ 
tached to the reflecting surface, the (local) incident frame 
K, = (0,Xi,y, Zi) attached to the incident beam, and 
the (local) reflected frame K r = (0,x r ,y, z r ) attached 
to the reflected beam. These three reference frames are 
connected via a rotation of an angle 6 around the y di¬ 
rection [37]. The reflecting surface is located at z = 0 , 
with the z-axis pointing towards the interface. With this 
choice of geometry, the incident beam comes from the 
region z < 0 and propagates in the x-z plane. 

The electric field in the incident frame can be then 
written, using its angular spectrum representation [36], 
as follows: 


E i( r ) = X f d 2 Ke x (U,V,d)A x (U,V,d)e ik ' r ', (2) 

A=1 ^ 


where d 2 K = dUdV, e x (U,V,6) is the local reference 
frame attached to the incident held [38], A X (U,V,9 ) = 
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Fig. 1. (Color online) (a) Schematic representation of the 
considered surface. The graphene layer (red) is charac¬ 
terized by its surface conductivity er(ft), whose explicit 
expression is given by Eq. (1). The dielectric substrate 
(green) is characterized by the refractive index n. (b) 
Geometry of beam reflection at the interface. The sin¬ 
gle graphene layer is located on the surface at z = 0. 
The different Cartesian coordinate systems K, K,. K r are 
shown. 


a\(U, V, 0)A{U, V ) and k, • r, = UX t + VY t + WZ U be¬ 
ing Xi = koXi the normalized coordinate in the inci¬ 
dent frame. Y) and Zi are defined in a similar manner. 
a\{U. V. 0) = e\(U,V,0) ■ f accounts for the projection 
of the beam’s polarization f = / p x + f s y (normalized 
according to |/ p | 2 + |/ s | 2 = 1) onto the local basis, and 
A(U, V) is the beam’s spectral amplitude, which here is 
assumed to be Gaussian, i.e., 


cos 0 — \J n 2 — sin 2 9 — er(ft) 

cos 9 + \/ n 2 — sin 2 9 + er(ft) 

n 2 cos 9 — \Jn 2 — sin 2 9 [1 — cr(ft) cos 9] 

n 2 cos 9 + \/n 2 — sin 2 9 [1 + cr(ft) cos 9] 


(5a) 


(5b) 


where 9 is the incident angle, n is the refractive index 
of the dielectric medium and er(ft) is the graphene’s sur¬ 
face conductivity, as defined by Eq. (1). The modulus R\ 
and phase <j>\ of the reflection coefficients r\ = R\e l< ^ x 
(with A £ { p , a}) are shown in Fig. 2, together with 
the correspondent quantities for the case of a simple di¬ 
electric surface without the graphene coating. While the 
presence of the graphene layer does not modify significa- 
tively R\ for neither p- or s-polarization (as it appears 
clear from Figs. 2(a) and (c), respectively), the change 
induced in the phases <j>\ of the reflection coefficients is 
considerable. For a normal air-glass interface, in fact, we 
have d<f>x/d9 = 0 being 9 the angle of incidence. Here, 
instead, we have 8cj)\/d9 ^ 0. A closer inspection of Eqs. 
(5), moreover, reveals that such a novel 0-dependence of 
the phases (f>\ is entirely due to the graphene conductiv¬ 
ity <7 (ft). 

To compute the GH and IF shifts, we calculate the cen¬ 
ter of mass of the intensity distribution in the reflected 
frame, namely [38] 


(R> =' 
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= (X r )± r + (Y r ) yy, (6) 


where R = (X r ,Y r ) T - Spatial (A) and angular (0) GH 
and IF shifts are then defined as follows: 


A(U,V) = e -^ (c/2+y2) , (3) 

being Wq the spot size of the beam. In the remaining 
of the manuscript, we will consider only well collimated 
beams, namely the paraxial assumption U, V <C 1 is im¬ 
plicitly understood. 

Upon reflection, the electric field can be then written 
as follows: 


E r (r r ) = ^ f d 2 Ke\(—U,V, tt — 0)A\(U,V, 0)e ikr ' rr , 

A=1 ^ 

( 4 ) 

where k r • ry = —UX r + VY r + WZ r and A\(U , V, 9) = 
r\(U, V, 9)A\(U , V, 9), with r\(U, V, 9) being the Fresnel 
reflection coefficients associated to the single plane wave 
component of the field [39]. The minus sign in front of 
U in e\, as well as in k r • ry accounts for the specular 
reflection of the single plane wave component [38]. 

The presence of a single layer of graphene deposited on 
the dielectric surface modifies its reflection coefficients as 


A gh = (X r ) , 

z=0 

A if = (Y r ) , 
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The explicit expressions of the GH and IF shifts for a 
fundamental Gaussian beam read, according to [41], as 
follows: 
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where f p = a p , f s = a s exp (ip) and w\ = a\R^/(a 2 R 2 + 
a 2 R 2 ) (where A £ {p, s}) is the fractional energy con¬ 
tained in each polarization state. 
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Fig. 2. Modulus (left column) and phase (right column) of the reflection coefficients r\ = R\ exp for p- 

polarization (top row) and s-polarization (bottom row). In all graphs, the solid black line corresponds to the case 
of the graphene-coated surface, while the red dashed curve corresponds to the case without graphene coating. The 
refractive index of the bulk medium is chosen to be n = 1.5. 


As suggested by Figs. 2(a) and (c), the changes in 
R\ introduced by the graphene layer are negligible. We 
therefore expect to observe no changes in the angular 
shifts €>gh and 0/f, as they are functions of R\ solely. 
The spatial shifts A gh and Ajp, on the other hand, 
contain a dependence on the phases (j )\, and they are 
therefore affected by the presence of the graphene coat¬ 
ing. Let us first discuss the IF shift. In this case 4> p — cf> s 
is very close to tt [Figs. 2(b) and (d)], and the resulting 
spatial shift A jp will be nonzero (but very small) even 
for linear polarization, in contrast to the case without 
graphene. 

More interesting is the case of the spatial GH shift. 
For a normal air-dielectric interface, one has d(f>\/dd = 0 
and therefore, according to Eq. ( 8 a), A gh = 0. It is in 
fact well known since the pioneering work of Goos and 
Hanchen [1], that A gh 0 occurs only in total internal 
reflection, where R\ = 1 and d(f)\/d9 ^ 0. For the case of 
a graphene-coated surface, on the other hand, the phase 
cf>\ varies with 9 for both s- and p-polarizations, as Figs. 
2 (b) and (d), respectively, show. In this case, then, we 
observe a nonzero spatial GH shift even without total 
internal reflection. 

The spatial GH shift A gh occurring at a graphene- 


coated dielectric surface is depicted in Fig. 3(a) and (b) 
for p- and s-polarization, respectively. As can be seen, 
for both polarizations we have A gh 7 ^ 0 although no to¬ 
tal internal reflection takes place. In particular, (f> p varies 
very rapidly from 0 to — 7 r in the vicinity of the Brewster 
angle 9b- This corresponds to a giant and negative spa¬ 
tial GH shift. On the other hand, <p s varies very smoothly 
with 9, thus resulting in a nonzero (but very small) spa¬ 
tial GH shift for s-polarization. 

In conclusion, we have presented a detailed theoretical 
analysis of GH and IF shifts of a Gaussian beam imping¬ 
ing onto a graphene-coated dielectric surface. Our anal¬ 
ysis revealed that the main effect of the graphene layer 
is to introduce, through its surface conducibility 
a dependence of the phases (j)\ of the reflection coeffi¬ 
cients on the incident angle 9. This, ultimately, reflects 
in the appearance of a nonzero spatial GH and IF shifts. 
In particular a giant and negative spatial GH shift in 
the vicinity of the Brewster’s angle for p-polarization has 
been predicted, in ageeement with the recently published 
experimental results [33]. The authors thank the German 
Ministry of Education and Science (ZIK 03Z1HN31) for 
financial support. 
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Fig. 3. Spatial GH shift A gh for (a) p- polarization and (b) s-polarization for a graphene-coated surface. Since 
d(f>\/d0 ^ 0, in both cases A qh ^ 0. In particular, since <f> p varies very rapidly with 9 in the vicinity of the Brewster 
angle, the corresponding spatial shift for p-polarization [Panel (a)] is giant in modulus, and negative due to the fact 
that (j) p varies from 0 to — 7r [See Fig. 2(b)]. 
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